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ABSTRACT 

The  equation  of  motion  for  an  elastically  mounted,  linearly  damped 
circular  cylinder  subjected  to  a  time-dependent  transverse  fluid  force 
has  been  formulated.  The  steady-state  and  transient  solutions  have 
been  obtained  in  terms  of  the  in-phase  and  out-of-phase  force  coeffi- 
cients, amplitude  ratio,  and  the  internal  damping  of  the  cylinder.   The 
results  have  shown  that  a  cylinder  may  be  self-excited  in  the  range  of 
Strouha I  frequencies  close  to  the  natural  frequency  of  the  cylinder. 
It  has  further  been  shown  that  the  heuristic  models  previously  used  are 
inadequate  to  correctly  predict  the  in-phase  or  the  inertia  component 
of  the  exciting  force.  Additional  experiments  and  analysis  could  en- 
able one  to  predict  the  behavior  of  cables  subjected  to  ocean  currents. 
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I.   INTRODUCTION 

Experience  has  shown  that  elastic  structures  of  one  or  more  degrees 
of  freedom,  such  as  tall  slender  structures,  cables,  etc.,  can  extract 
energy  from  the  flow  about  them  and  can  develop  catastrophic  flow-induced 
oscillations.   The  understanding  of  this  energy-extraction  process  is  of 
paramount  importance  if  one  is  either  to  eliminate  or  minimize  it  or  to 
design  the  elastic  structure  such  that  it  can  withstand  the  oscillations 
under  the  contemplated  environmental  conditions. 

Evidently,  the  real  origin  of  the  self-excited  oscillations  lies  in 
the  basic  mechanism  of  the  vortex-shedding  process.  Consequently,  the 
prediction  of  such  oscillations  requires  a  tractable  theoret ica I -formu I a- 
tion  of  the  vortex-shedding  phenomenon.  The  fact  that  such  a  formulation 
does  not  exist  even  for  separated  flow  about  a  stationary  cylinder  shows 
the  complexity  of  the  treatment  of  flow  about  oscillating  bodies.  The 
problem  would  have  been  relatively  simpler  had  the  vortex  shedding  and 
the  forces  arising  from  it  been  independent  of  the  motion  of  the  body. 
Experiments  show  that  vortex  generated  forces  give  rise  to  oscillations, 
and  the  oscillations  modify  the  forces,  and  so  on,  finally  resulting  in 
a  steady-state  oscillation  such  that  the  Strouha I  frequency  of  the  vortex 
shedding  nearly  coincides  with  the  natural  frequency  of  oscillation  of 
the  body. 

Mathematical  models  of  flow- induced  vibrations  of  bluff  bodies  and 
the  response  of  circular  cylinders  to  vortex  shedding  have  been  aptly 
described  by  Parkinson  [Jj,  and  Currie,  et  al.  [2~]   and  will  not  be 
repeated  here. 


It  appears  that  among  the  various  models  considered  so  far,  the 
"wake-oscillator"  model  of  Hartlen  and  Currie  [3]]  attracted  considerable 
attention  among  the  students  of  vibration  analysis.   This  model  predicts 
some  of  the  observed  features  of  the  self-excited  oscillations  even 
though  the  reasons  as  to  why  an  osci  I lator  of  the  van  der  Pol  type 
should  at  all  describe  the  self-excited  oscillations  of  an  elastic  body 
are  not  yet  clear.  There  are  several  coefficients  in  the  model  which 
must  be  determined  experimentally.   There  are,  at  present,  wide  varia- 
tions in  the  methods  of  measurement  of  some  of  the  coefficients.   For  ex- 
ample, the  damping  of  the  vibrations  of  a  given  beam  or  cylinder  is 
determined  in  still  air  or  water  depending  on  whether  the  self-excited 
oscillations  take  place  in  air  or  water.  The  added  mass  is  assumed  to 
be  equal  to  the  displaced  mass,  of  the  cylinder.  There  is  no  distinction 
between  the  material  damping  as  would  be  determined  in  vacuum,  and  the 
damping  caused  by  the  fluid. 

It  became  clear,  during  the  course  of  the  present  investigation, 
that  the  forces  acting  on  a  cylinder  oscillating  harmonically  transverse 
to  a  uniform  flow  should  be  measured  and  the  in-phase  and  out-of-phase 
components  of  the  force  should  be  determined.  This  phase  of  the  inves- 
tigation has  been  carried  by  Meyers  [4]  and  Fortik  \J5~]   under  the 
direction  of  Professor  T.  Sarpkaya,  and  the  drag  and  inertia  coeffi- 
cients for  the  lift  force  have  been  determined  through  the  use  of 
Morison's  [_6^\   equation. 

The  purpose  of  the  work  reported  herein  is  to  make  use  of  the 
measured  force  by  combining  it  with  the  equation  of  motion  of  an 
elastically  mounted,  linearly  damped,  circular  cylinder  and  to  attempt 
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to  predict  the  self-excited  oscillations  of  the  cylinder  for  given  and 
directly  measurable  physical  parameters  such  as  the  mass,  length, 
diameter,  material  damping,  and  the  natural  frequency  of  the  body,  and 
the  density  and  velocity  of  the  fluid. 

It  will  be  assumed  that  the  measurements  and  the  predictions  are 
in  the  same  Reynolds  number  range  and  that  the  force  coefficients  do 
not  measurably  depend  on  the  Reynolds  number  in  the  range  under  considera- 
tion. 
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I  I.   ANALYSIS 

A.   FORCING  FUNCTION 

The  transverse  force  acting  on  a  circular  cylinder  undergoing  har- 
monic oscillations  in  the  transverse  direction  to  a  uniform  flow  is 
expressed  in  terms  of  the  Morison's  equation  as  C-O 

c  =  ,   F  ,  =  c  .  tt2  (JhI)  A2  sin  *T  t 
L   iPfLDV2    ml      D    VT       T 

_   ,  m  ,2  ,D,2  i    2tt  ,  I     2tt  ,  ,  . 

-  C  . .  (— k-)   ( — )   cos  -=-  t  cos  -=-  t  (  I  ) 

d  I   U     ^T   '     I   '      I 

where  C  ,  and  C,.  are  the  inertia  and  drag  coefficients,  respectively, 
ml      d  I 

Evidently,  the  normalized  force  as  well  as  the  drag  and  inertia  co- 
efficients depend  on  the  parameters  U  T/D  =  2ttA/D,  and  D/VT.   The  first 

m 

parameter  is  directly  proportional  to  the  relative  amplitude,  and  the 
second  parameter  may  be  regarded  as  the  frequency  parameter  (D/VT)  or 
the  reduced  velocity  (VT/D). 

For  the  purpose  under  consideration  it  is  more  advantageous  to  ex- 
pand the  |  cos  cot  |  cos  cot  term  in  series  and  to  retain  only  the  first 
term.  This  procedure  yields, 


C.    =  C    .    it    (— =— )  ( — )      sin  cot  -  ■=—  C  , ,  (-7T-)    ( — )      cos  cot  (2) 

L         m  I  D       ttt  3tt     d        D  yT 


Equation  (2)  could  be  simplified  further  by  combining  it  with  the  equa- 
tion of  motion  for  an  elastically  mounted  cylinder. 
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B.   ELASTICALLY  MOUNTED  CYLINDER 


Fig.  I.  Elastically  mounted,  linearly 
damped  cy I  i  nder. 


The  equation  of  motion  for  the  elastically  mounted,  linearly  damped, 
and  periodically  forced  cylinder  may  be  written  as  (see  Fig.  I) 


mx  +  ex  +  kx  =  F  =  j   p  DL  V  C 


(3) 


where  m  represents  the  mass  of  the  cylinder,  c  the  linear  damping,  and 

k  the  spring  constant.  The  derivatives  of  x  are   taken  with  respect  to 

...      .     d  x 
real  time  t,  i.e.,  — *■  =  x,  etc. 

dt 

Def in  ing, 


oj  =  /k/m  ,  oj 

m       '   n 


2tt  f   ,  t  =  oj  t  ,  x/D 
n        n 


=  x 


5  = 


2m  oj  '  o 


n  =  f  /f  ,  a 

c  n  ' 


f  /f   =  SV/f  D  ,   X 
s  n       n      r 

2 
.  LDZ  f 
pf      s 


2  2  2' 
2m  oo  S  oj 
n    o 


=  A/D 
3  x 


3t' 
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equation    (3)   may   be   reduced   to 

?  ?  OS   2 

x     +2#     +  x     =   aw  [C    ,    7T   (2ttX    )  (— )      sin   fix 
r  r         r  o     ml  rco 

o 


^-C^,(2ttX    )2    (— )2   cos   Ox]  (4) 

3tt     d  I  r  co 

o 


3      2 
Defining   p     =    p,/p     and   noting   that   p     =   2tt  aS    ,    equation    (4)    re- 


duces to 


"?  1  f~\ 

x     +  2£<     +  x     =   p  X  ft   (C    .    sin   fix *■  X     C. .    cos   0T)  (5) 

r  r  r         r   r         ml  _  2     r     dl 

3tt 


or   by   a   change  of   variable   from  x     to  y     =  x/A,    and    further   simplifica- 
tion one   has 

y     +  2Cy     +  y     =   P  ft2(C    ,    sin   Qt  -  -^4-  X     C  . .    cos   fix)  (6) 

7  r  r        '  r  r  ml  7   2      r     d  I 

3tt 


Equation  (6)  may  be  written  as 

y  +  2£y  +  y  =  EQ,     sin  (Ox  -  a)  (7) 

r      r    r 


in  which 


/ 


E  "  Pr  V  Cm|2  +  <^>(B»  Cd|:2  C8J 


and 

I6A  C 
a  =  arc  tan  * (9) 

3ir  D  C  . 

ml 

2 
The  coefficient  EQ,     may  also  be  expressed  in  terms  of  the  maximum 

lift  coefficient  defined  by 

F  =  I  p  V2DL  (C, )   sin  (0T  -  a)  (10) 

Z      t  L  m 
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or 


Inserting  equation  (10)  into  equation  (3)  and  simplifying  one  has 

I    -2 
n       4r   px  V  DL(C,  ) 

EP.2  =  111 _ LiH  en, 

m  A  oj 

n 


VT  2 
C-tt)    (Ci  } 

E  =  pr— T-  WDT  (l2) 

2tt 


Combining  equations  (12)  and  (8)  and  simplifying  one  has 


CL}m  ,   2tt3    /  "  2     16   k  r   .2  fm 

D 


It  should  be  noted  that  if  (A/D)  C  .  were  independent  of  A/D,  then 

it  is  easy  to  show  that  (C.  )  /(A/D)  would  also  be  independent  of  A/D 
'  L  m 

since  the  experiments  have  shown  that  C   is  independent  of  A/D.   This 

ml 

would  make  the  results  independent  of  A/D,  a  conclusion  which  is  contrary 

to  observations.  This  proves  that  (A/D)  C_,,  and  (C,  )  /(A/D)  are  indeed 

r  d  I       L  m 

dependent  on  A/D  in  a  complex  manner.   But  this  dependence  is  rather 
weak.  That  is  why  slight  changes  in,  for  instance,  p  can  lead  to  large 
changes  in  £,  or  Q. 


C.   SOLUTION  OF  THE  EQUATION  OF  MOTION 

As  shown  above,  a  convenient  form  of  the  equation  of  motion  is 

y  +  2£y  +  y  =  Eft2  s i n  (Or  -  a)  (14) 

7  r  r        7  r 

I  .  Transient  Solution 

The  Homogeneous  equation  corresponding  to  equation  (14)  is 

V\  +  2£yr  +  y  =  0  (15) 

r      r    r 
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Its  characteristic  equation  reduces  to 


s  +  2?s  +  I  =  0 


(16) 


whose  roots  are  given  by 


,2 


=  -c  ±  i/l  -  C: 


(17) 


Considering  only  the  underdamped  case,  it  follows  that  £  <  I. 
Therefore,  according  to  Boyce  and  Di  Prima  C7~J,  the  solution  to  the 
corresponding  homogeneous  equation,  that  is,  the  transient  solution  is 
given  by 


(y  ),  =  Y,  e~CT  sin  (/l  -  C2  T  -  8) 
7r  h    h 


(18) 


where  Y,  and  3  will  be  determined  from  the  initial  conditions, 
h 

20  Steady-State  Solution 

The  particular  solution  to  equation  (14)  is  a  steady-state 

oscillation  of  the  same  frequency  of  the  excitation.  Thus,  as  indicated 

in  Thomson  \j3i2t    i"he  particular  solution  is  assumed  to  be  of  the  form 


(y  )  =  Y  sin  [(flr  -  a)  -  <j>] 
r   p    p 


(19) 


where  (J)  is  the  phase  of  the  displacement  with  respect  to  the  excitaTion, 
as  shown  in  Fig.  2,  and  is  given  by 


Fig.  2.  Vector  diagram. 
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and 


$  =  arc  tan  2#  (20) 

I -ft2 


and  Y   is  qiven  by 

P 


Eft2 


Y  =  =         ■  (21) 

p   ±  /(I -ft2)2  +  (2<;ft)2 


Adding  the  transient  and  the  steady-state  solutions,  the  following  com- 
plete solution  is  obtained 


y  =  Yu  e"^T  sin  (/l-<;2  T+g)  +  Y  sin  (flx-a-<fr)  (22) 

rr    h  p 


Also,  its  derivative  with  respect  to  t  is  given  by 


y     =  Yu{-^e"CT   sin    (/l-£2   x+3)    +  e   CT   /l-?2   t+6)} 
r  h 

+  Y     ft  cos(ftT-cH>)  (23) 

P 

Applying  the  first  initial  condition  that  x  =  0  corresponds  to  y  =0, 
one  finds 

Y  sin  C-a-4)  (24) 

Y  =  -  -E — 

h        sin  6 

Applying  the  second  initial  condition  that  t  =  0  corresponds  to  y  =  0, 
one  obtains 


0  =  Y.C-5  sin  6  +  /l-^2  cos  6)  +  Y  ft  cos(-OMfr)   (25) 
nserting  the  expression  for  Y  given  by  equation  (24)  in  equation  (25), 
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one  has 


.2 


6  =  arc  tan  (   *'  %, jt-  (26) 

C+ft  cot(-a-<j>) 


A I  so,  reca  I  I  i  ng  that  as  t  ->  °°,  x  ■*■  -A  sin  $,  and  that  y  =  x/A  •>  -sin  9 , 
it  is  seen  from  equation  (22)  that  Y  =  -I.   Introducing  this  value  in 
equation  (21)  and  noting  that  E  >  0,  one  obtains 


=  Ai-ft2)2  +  (2ffi)2  27 

ft2 


Substituting  the  expressions  for  Y  and  Y   into  equation  (22),  one  ob- 
tains the  following  complete  solution  of  the  equation  of  motion 


sin(-a-4>)  -£t  .  ,  J,    ~7  ,ON   .  fr,  ,.  /on, 

y  =  : — tt— e         sin  (/ -£  x+8)  -sin  (flr-a-A)  (28) 

'  r     si  n  3  *  r 


The  results  obtained  through  the  use  of  equation  (28)  will  be 
presented  in  connection  with  the  discussion  of  results. 
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III.   DISCUSSION  OF  RESULTS 

The  time  dependent  force  acting  on  the  cylinder  was  expressed,  as 

indicated  previously,  in  terms  of  C  .  and  C...  Assuminq  this  force  to 
r        '  ml      d I 

be  applied  to  an  elastically  mounted  and  linearly  damped  cylinder,  the 

equation  of  motion  for  the  cylinder  was  developed  and  solved  in  terms 

of  the  given  parameters  and  experimentally  obtained  coefficients.  Two 

of  these  coefficients,  C  .  and  C_,.,  inertia  and  draq  coefficients  respec- 

ml      d  I 

tively,  were  introduced  in  the  expression  of  the  time  dependent  force 

acting  on  the  cylinder.   Using  the  experimental  apparatus  and  procedure 

described  in  C4H,  data  were  obtained,  which  together  with  the  results 

from  r4l.  were  used  to  plot  C  .  and  C . .  as  a  function  of  VT/D  for 

r     ml      d  I 

A/D  =  0.25,  0.50,  0.75  and  0.84,  as  shown  in  Figs.  3  and  4. 

The  results  plotted  in  Fig.  3  show  that  for  the  region  of  interest, 

i.e.,  for  VT/D  in  the  range  of  4  to  7,  and  for  the  precision  of  the 

experiments,  C  ,  is  independent  of  A/D. 
ml 

On  the  other  hand,  C  .  is  dependent  on  both  VT/D  and  A/D  as  shown 

in  Fig.  4.  As  noted  previously,  (A/D)  C_, ,  and  (C.  )  /(A/D)  are  not  con- 
3  "        d  I       L  m 

stants.   Furthermore,  it  can  be  shown  that  (C,  )  /(A/D)  decreases  as 

L  m 

A/D  increases.,   Figures  5,  6,  7  and  8  show  (C,  )   as  a  function  of  VT/D 

3      '   '  L  m 

for  A/D  =  0.25,  0.50,  0.75  and  0.84  respectively. 

Therefore,  an  iterative  procedure  is  needed  to  determine  the  complete 

solution  of  the  equation  of  motion  for  a  specific  case,  as  indicated 

below. 

Firstly,  for  a  given  VT/D,  C  ,  is  obtained  from  Fiq.  3.  Then,  for 
i  r  z>  'ml  3  ' 

a  specific  relative  amplitude  A/D  and  the  same  VT/D,  C  .  is  obtained  from 
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Fiq.  4,  and  the  value  of  (C.  )   is  obtained  from  Figs.  5,  6,  7  or  8,  or 
3  L  m 

alternatively  from  equation  (13).  Then,  E  is  calculated  using  equa- 
tion (12),  since  p   is  known.  Substituting  this  value  for  E  in  equa- 

^Tn   VT 
tion  (27),  a  value  for  £  is  calculated,  since  ft  =  (— =— )/(-?r-) .  A  plot 

of  l/E  versus  ft  for  specific  values  of  £  is  shown  in  Fig.  9.  Then, 

the  values  for  a,  $  and  3  are  calculated  from  equations  (9),  (20)  and 

(26),  respectively.  Substituting  these  results  into  the  equation  of 

motion  and  plotting  the  results,  one  obtains  curves  similar  to  those 

shown  in  Figs.  10,  II  and  12.  These  figures  show  clearly  the  effect  of 

£  on  the  response  of  the  system.  A  lowly  damped  system  results  in  an 

initial  overshoot  of  the  amplitude  before  it  reaches  a  steady-state. 

A  highly  damped  system,  on  the  other  hand,  reaches  its  steady-state 

gradually  from  its  initial  rest  position. 

Evidently,  one  can  perform  various  parametric  studies  on  the  effect 
of  damping,  density  ratio,  etc.,  on  the  transient  and  steady-states  of 
the  oscillation  provided  that  sufficient  experimental  data  are  avail- 
able to  obtain  accurately  enough  the  drag  and  inertia  coefficients.   It 
is  also  evident  that  the  practical  problem  is  the  determination  of 
whether  a  cylinder  of  given  natural  frequency,  density  ratio,  internal 
damping,  etc.,  will  be  self  excited  in  a  fluid  of  given  velocity. 

The  determination  of  this  problem  requires  an  iteration  not  only 
over  VT/D  but  also  over  A/D  in  order  to  obtain  the  appropriate  values 
of  the  drag  and  inertia  coefficients.  Preliminary  attempts  to  resolve 
this  aspect  of  the  problem  have  resulted  in  the  computer  program  given 
in  Appendix  A.  This  program  needs  to  be  improved  with  additional  data 
particularly  to  the  lower  and  higher  values  of  A/D  where  the  force 
transfer  coefficients  vary  most  rapidly  with  A/D. 
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Then,  it  will  be  possible  to  improve  the  program  to  predict  with 
greater  precision  the  range  of  self-excited  oscillations  as  well  as 
their  amp  I i  tude. 

A  computer  program  for  the  complete  solution  given  by  equation  (28) 
and  for  the  plotting  of  the  results  is  given  in  Appendix  B. 

It  is  evident  from  the  foregoing  discussion  that  the  prediction  of 
the  characteristics  of  the  self-excited  oscillations  is  based  on  the 
measured  time-dependent  transverse  force  rather  than  on  an  assumed 
variation  of  the  lift  coefficient  through  the  use  of  a  heuristic  model 
C33.   It  was,  therefore,  desirable  to  compare  the  components  of  the 
alternating  force  obtained  from  the  experimental  data  with  those  pre- 
dicted by  the  Hartlen  and  Currie  [J>~]   model.   For  this  purpose  the  i  n- 
phase  and  out-of-phase  components  of  the  Hartlen  and  Currie  transverse 
force  coefficients  have  been  written  equal  to  the  corresponding  com- 
ponents of  the  model  used  herein  [see  equation  (2)3.  This  procedure 
y iel ded 

32tt  (^)2 

C  .(Hartlen  and  Currie)  = —  C,,  (29) 

d  VT  2   dl 

3(— ) 
D 

and 

C  (Hartlen  and  Currie)  = —  C  ,  (30) 

m  VT  2    ml 

(— ) 

D 

in  which  C  .  and  C  , .  are  the  drag  and  inertia  coefficients  evaluated  in 
ml      d  I 

the  present  investigation. 

Equations  (29)  and  (30)  have  been  evaluated  for  A/D  =  0.25  and  0.50 

and  the  results  are  plotted  as  a  function  of  f  /f   in  Fiqs.  13  and  14. 

c  s      3 
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A  similar  plot  obtained  from  [33  is  shown  in  Fig.  15  for  A/D  =  0.2. 

Even  though  the  A/D  values  are  not  exactly  identical,  Figs.  13,  14  and 

15  could  certainly  be  compared  as  far  as  the  overall  variation  of  C 

and  Cj  are  concerned.   It  is  evident  that  C,  curves  are  fairly  similar 
d  d 

and  that  the  Hartlen  and  Currie  model  could  be  adjusted  sufficiently  to 

yield  similar  C,  values  for  identical  A/D  values.  However,  C  values 
;  d  '  m 

which  are  comparable  only  for  f  /f  <  I  show  opposite  variations.   It 
has  been  previously  stated  that  the  use  of  an  added  mass  coefficient 
in  terms  of  the  displaced  mass  of  the  body  is  not  sufficient  and  that 
the  inertia  coefficient  of  a  body  oscillating  in  the  transverse  direc- 
tion to  a  uniform  flow  could  be  significantly  different  from  its  dis- 
placed mass.   It  suffices  to  conclude  that  the  Hartlen  and  Currie  model 

does  not  predict  C   in  accordance  with  the  observations  and  that  the 
r       m 

overall  predictions  of  the  model  may  be  attributed  partly  to  the  adjust- 
ment of  the  other  free  parameters  such  as  the  damping  ratio,  the  density 
ratio,  etc. 
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IV.   CONCLUSIONS 

The  results  presented  herein  have  shown  that  (a)  a  cylinder  can 
undergo  self-excited  oscillations  for  certain  values  of  the  ratio  of 
the  natural  frequency  of  the  cylinder  to  the  Strouha I  frequency  of  the 
vortices  (f  /f  ),  internal  damping  of  the  cylinder  (£),  and  the  density 
ratio  (p  );  (b)  the  drag  and  inertia  coefficients  must  be  determined 
experimentally  as  a  function  of  A/D  and  VT/D  and  must  be  incorporated 
into  the  equation  of  motion;  (c)  the  prediction  of  the  self-excited 
oscillations  and  their  amplitude  will  be  possible  with  the  acquisition 
of  additional  data  for  the  transverse  force;  and  that  (d)  the  results 
presented  herein  are  quite  encouraging  and  certainly  a  reflection  of 
the  combination  of  the  fluid  dynamical  and  vibrational  aspects  of  the 
phenomenon  unlike  the  heuristic  models  previously  proposed. 
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Scale: 
J2  :   0.6/in 
l/E:  0.5/in 


Fig.  9.  The  function  l/E  versus  Q  for  various  values  of  the 
damping  ratio  5. 
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Fig.  II.  Sample  calculations  for  the  transient  and  steady-state 
oscillations  of  a  cylinder  for  the  specific  parameters 
ind  icated. 
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Fig.  12.  Sample  calculations  for  the  transient  and  steady-state 
oscillations  of  a  cylinder  for  the  specific  parameters 
i  nd  icated. 
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Fig.  15.  C  and  C  versus  f  /f  for  A/D  =  0.2,  from 
Hartlen  and  Currie    C3H. 
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